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EXECUTIVE  SUMMARY 


Field  studies  were  conducted  on  the  Holston  River  in  the  vicinity 
of  Holston  Army  Ammunition  Plant  (HAAP)  near  Kingsport,  Tenn.  to  obtain 
hydrological  and  environmental  parameters  necessary  for  the  computer 
modeling  of  the  fate  (loss  and  movement)  of  the  munitions -related  water 
pollutants,  2,4,6-trinitrotoluene  (TNT)  and  hexahydro-l,3,5-trinitro- 
1,3,5-triazine  (RBX) .  Holston  River  water  samples  were  collected  at  0, 

4.5,  7.9  and  33.0  km  from  the  discharge  point  and  analyzed  for  TOT  and 
RDX  as  a  means  of  validating  the  concentration-time  dependent  predictions 
of  the  computer  model. 

A  one-dimens-*  *>nal  compartmentaJ  Ized  river  system  was  found  to 
accurately  predict  RDX  concentrations  at  distances  greater  than  20  km 
from  the  discharge  point  where  the  HAAP  discharge  stream  becomes  uniformly 
distributed  in  the  Holston  River.  Up  to  20  km,  the  waste  stream  is  not 
uniformly  distributed  and  a  two-dimensional  compartmentalized  river 
system  was  developed  to  account  for  the  lateral  dispersion  effect  of 
HAAP  waste  streams  in  the  Holston  River.  The  transformation  of  RDX  by 
environmental  processes  was  determined  not  to  be  significant  over  the 
studied  river  reach  (33  km),  and  therefore,  the  two-dimensional  model 
network  describes  the  dilution  processes  in  the  Holston  River. 

A  pseudo-first-order  photolysis  rate  constant  of  4,0  day  ^  and 
pseudo-first-order  biotransformation  rate  constant  of  0,03  day  ^  for 
TNT  were  estimated  from  laboratory  studies  and  field  data.  Using  the 
two-dimensional  model,  concentration-time  dependent  curves  were  constructed 
to  predict  the  concentration  of  TNT  as  a  function  of  distance  and  location 
(North  Bank,  Center  of  River,  and  South  Bank)  in  the  Holston  River.  Due 
to  the  low  input  of  TNT  into  the  Holston  River  during  the  field  study, 
sufficient  data  could  not  be  collected  to  validate  the  model  for  TNT. 

We  were  able,  however,  to  determine  that  the  laboratory  estimated  photo- 
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chemical  rate  constant  was  consistent  with  that  measured  in  the  Holston 
River. 

Field  study  data  collected  for  the  compounds  1-acetyloctahydro- 
3,5,7-trinitro-l»3,5 ,7-tetrazocine  (SEX) ,  l-acetylhexahydro-3,5-dinitro 
1,3,5-triazine  (TAX),  and  octahydro-1,3,5 ,7-tetranitro-l,3,5 ,7-tetra~ 
zocine  (HMX) ,  suggest  that  these  munitions-related  water  pollutants  do 
not  transform  readily  in  the  Holston  River. 
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I  INTRODUCTION 


the  U.S.  Army  is  in  the  process  of  assessing  the  hazards  associated 
with  muni t ions -related  water  pollutants  in  the  aquatic  environment .  Fart 
of  this  assessment  is  developing  an  understanding  of  the  lass  and  movement 
(fata)  of  pollutants  when  acted  upon  by  various  environmental  processes. 

In  Phase  II  of  our  environmental  fate  assessment  program  (Spanggord  et 
al. ,  1980b),  we  identified  the  dominant  transport  and  transformation 
processes  affecting  the  persistence  of  2,4,6— trinitrotoluene  (TNT),  hexa- 
hydro— 1,3, S-trinitro— 1 ,3, 5-triazine  (RBX),  2?4~diBitrotoluene,  and 
trinitrcglycerin  la  the  aquatic  environment.  We  also  applied  the  SRI 
aquatic  computer  model  to  selected  water  bodies  receiving  munition  wastes 
to  generate  concentrations  of  pollutants  as  a  function  of  time,  which 
would  allow  the  estimation  of  pollutant  concentrations  at  various 
locations  downstream  from  the  discharge  point  into  these  water  bodies. 

The  SSI  International  aquatic  fate  model  can  be  used  to  estimate 
pollutant  concentrations  reaching  downstream  water  intake  supplies,  to 
estimate  exposure  levels  to  various  aquatic  and  mammalian  populations,  and 
to  estimate  pollutant  concentrations  at  plant  boundaries  in  the  event  that 
regulatory  agencies  establish  guidelines  for  the  pollutants  being  dis¬ 
charged.  Because  of  these  potential  uses,  it  becomes  important  to  esta¬ 
blish  the  reliability  of  the  model  to  predict  pollutant  concentrations 
under  a  given  mt  of  environmental  conditions. 

The  objective  of  this  study  was  to  validate  the  SRI  computer  model 
used  in  Phase  II  by  selecting  a  water  body  receiving  munitions  wastes, 
determine  the  various  environmental  input  parameters  from  field  measure¬ 
ments,  predict  pollutant  concentrations  at  selected  locations  in  the  water 
body,  and  compare  the  predicted  values  with  actual  values  obtained  during 
the  field  study.  The  chemicals  selected  for  investigation  were  TNT  and 


The  results  of  this  study  can  be  used  to  establish  the  degree  of 
confidence  in  the  SRI  or  similar  EXAMS  model  for  estimating  pollutant 
concentrations  and  the  extent  to  which  these  models  may  contribute  to  the 
overall  hazard  assessment  of  munition-related  water  pollutants. 


II  BACKGROUND 


A.  SRI  International  Model 

The  study  reported  herein  is  based  on  the  application  of  the  SRI 
International  Computer  Model  (hereafter  called  the  SRI  model)  to  estimate 
the  concentration  of  pollutants  in  the  Holston  River  emanating  from  the 
Hols ton  Army  Ammunition  Plant  (HAAP)  near  Kingsport,  Tennessee. 

The  basic  model  was  developed  under  contract  with  the  U.S.  Envi¬ 
ronmental  Protection  Agency  (Smith  et  al.»  1977)  and  was  applied  to  the 
Holston  River  with  some  modifications.  The  model  requires  information  on 
the  hydrologic  parameters  of  the  site  under  investigation  and  the  rate  and 
equilibrium  constants  for  the  major  transport  and  transformation  processes 
that  act  on  the  chemicals  of  interest.  The  hydrologic  parameters  are  de¬ 
termined  from  field  study  measurements,  and  the  rate  constants  are 
obtained  from  laboratory  studies  using  the  natural  water  from  the  site. 

The  simplest  form  of  the  SRI  model  used  to  describe  a  flowing  river 
is  the  one-dimensional  compartmentalized  river  system.  A  diagram  of  this 
version  of  the  model  is  shown  in  Figure  i.  In  this  case,  a  section  of  the 
river  is  described  as  one  compartment  with  fixed  dimensions  of  length, 
width,  and  depth.  This  compartment  (water  compartment)  interacts  with  one 
sediment  compartment  directly  below  it.  The  chemical  disappears  at  an 
overall  rate  controlled  by  the  major  kinetic  and  equilibrium  processes 
established  for  a  specific  water  body.  A  river  simulation  Is  constructed 
by  joining  compartments  together,  and  the  loss  of  chemical  In  each  com¬ 
partment  is  determined  by  the  length  of  time  the  chemical  resides  In  the 
compartment  and  the  rate  constant  for  the  dominant  fate  processes.  The 
residence  time  is  calculated  from  the  river  flow  velocity.  Therefore,  the 
input  parameters  that  are  required  to  operate  the  model,  once  the  rate  and 
equilibrium  constants  have  been  measured,  include  length,  width,  and  depth 
of  the  compartment,  river  flow  velocity,  initial  chemical  loading,  and  pa¬ 
rameters  that  are  unique  to  specific  transformation  or  transport  processes 
such  as  light  intensity  (photolysis  and  oxidation),  microbial  population 
(biotransformation),  pH  (hydrolysis),  and  wind  speed  (volatilization). 


FIGURE  1  PHYSICAL  CONFIGURATIONS  OF  THE  COMPARTMENT  AND 
COMPARTMENTALIZED  RIVER  SIMULATIONS 


One  of  Che  major  assumptions  of  this  model  Is  that  each  compartment 
is  a  completely  aired  reactor.  To  account  for  cases  in  which  the  distri¬ 
bution  of  a  chemical  across  the  river  is  not  uniform,  we  employed  a  two- 
dimensional  compartmentalized  river  system  in  which  lateral  flow  veloc¬ 
ities  describe  the  movement  of  the  chemical  across  the  width  of  the  river 
(Figure  2).  In  this  case,  two  or  more  compartments  are  formed  across  the 
river  width.  The  movement  of  water  from  compartment  1  to  compartment  2  is 
equal  to  the  movement  of  water  from  compartment  2  to  compartment  1*  Lat¬ 
eral  flows  between  compartments  2  and  3  are  described  similarly* 
Therefore,  compartments  bordering  the  discharge  bank  will  dilute  con¬ 
tinually  until  homogeneity  is  achieved  across  the  river.  When  homogeneity 
is  observed,  through  either  dye  studies  or  the  monitoring  of  persistent 
chemicals  in  the  river,  then  the  lateral  flow  velocity  may  be  calculated 
from  distance  and  river  flow  velocity  measurements. 

B.  Mathematical  Description  of  Fate  Processes 

Based  on  our  Phase  IX  studies  (Spanggord  et  si.,  1980b),  the 
processes  that  should  be  considered  to  describe  the  loss  of  TNT  and  KDX 
are  photolysis,  biotransformation,  and  sorption. 

1.  Photochemical  Rate  Constant 

Photolysis  can  be  described  by  a  first  order  rate  expreasion 
(Equation  1) 


dt 


k  IC] 

pi  j 


(1) 


where  the  rate  of  loss  of  chemical,  C,  is  equal  to  a  pseudo-first-order 
rate  constant,  kp,  times  the  chemical  concentration,  [C]. 

The  rate  constant  is  a  function  of  the  sunlight  intensity,  which 
varies  with  the  time  of  day,  season,  latitude,  weather  conditions,  and  the 
light  scattering  properties  of  the  water  body.  To  account  for  the 
variations  in  kp,  the  rate  constant  measured  in  th?  laboratory  (sunlight) 
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at  a  particular  time  of  day  was  recalculated  by  a  model  program  (PHOTO) 
using  Equation  2 


k 


P 


K 

? 


sin  [ y-  (t-6)  |j  +  1.0 
- 5 - 


(2) 


where  kp  is  the  photolysis  rate  constant  at  noon,  and  t  is  the  time  of 
day.  This  equation  assumes  a  12-hour  sunlight  day  starting  at  6  AM  and 
ending  at  6  PM;  kp  reaches  a  maximum  at  12  noon  and  is  zero  between  6  PM 
and  6  AM.  A  plot  of  kp  versus  time  of  day  appears  in  Figure  3. 

2.  Biotransformation  Rate  Constant 

Biotransforaation  in  a  natural  water  body  is  believed  to  occur  under 
conditions  involving  a  large  microbial  population  and  a  low  chemical  con¬ 
centration,  Under  these  conditions,  the  rate  of  loss  of  chemical  can  be 
described  by  a  first-order  rate  expression  (Equation  3), 

~  ff  "  kb  <3) 


where  is  the  pseudo-first  order  rate  constant.  This  constant  is 
dependent  on  the  microbial  population.  To  account  for  cases  in  which  the 
microbial  population  varies,  the  model  employs  a  second-order  rate 
expression  shown  in  Equation  4. 


*Sr"kb2Wlcl 


(4) 


where  is  the  second-order  rate  constant  and  {X}  is  the  microbial 

population.  The  value  is  calculated  from  k^  and  [X]  using 
Equation  5. 
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3.  Sediment  Sorption  Partition  Coefficient 

The  sediment  sorption  partition  constant  is  obtained  from  the 
empirical  relationship  shown  in  Equation  6, 


C  -  kfC  ]l^° 
s  1  wJ 


(6) 


where  [C8]  is  the  amount  of  chemical  sorbed  per  mass  of  sediment,  Cw  is 
the  amount  of  chemical  per  unit  volume  of  water,  and  n  is  a  constant.  At 
equilibrium  and  low  chemical  concentrations,  n  is  often  nearly  equal  to 
1.  Under  these  conditions,  the  sediment  sorption  partition  coefficient 
is  defined  as  Equation  7. 


(7) 


K_  is  determined  in  the  laboratory  using  sediment  from  the  water  body 
6 

under  investigation  and  is  applied  directly  into  the  model. 

C.  Model  Input-Output  Computations 

Based  on  the  above  discussion,  the  model  input  parameters  Include 
compartment  volume,  the  photolysis  rate  constant  extrapolated  to  noon,  the 
second-order  biottansforaation  rate  constant  and  the  microbial  population, 
the  sediment  sorption  partition  coefficient,  and  the  chemical  loading  into 
the  first  compartment.  Sediment  settling  and  scouring  rates  (calculated 
from  sediment  loading  data  at  several  locations)  are  also  determined  to 
estimate  the  fraction  of  chemical  that  is  not  returned  to  the  aqueous 
phase • 
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From  these  data,  the  model  calculates  the  concentration  of  chemical 
transformed  and  the  distribution  of  cheaical  between  the  sediment  and 
aqueous  phases.  A  flow  chart  of  the  SHI  model  for  the  environmental 
assessment  appears  in  Figure  4. 

The  model  output  consists  of  changes  in  cheaical  concentration  and 
mass  in  the  aqueous  and  sediment  phases  of  each  compartment.  Cheaical 
concentration  at  any  point  from  the  discharge  point  in  the  water  body  can 
then  be  plotted  as  a  function  of  time. 
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FIGURE  4  FLOW  CHART  OF  THE  SRI  MODEL  FOR  ENVIRONMENTAL  ASSESSMENT 
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Ill  METHOD  0?  APPROACH 


A.  Site  Selection 

To  validate  the  SRI  model,  we  chose  to  study  the  Holstoa  River,  which 
receives  wastes  from  the  Kolston  Army  Ammunition  Plant  (HAAP)  in  Kings¬ 
port,  Tennessee*  This  river  was  selected  because  both  TNT  and  RDX  are 
discharged  simultaneously  into  it,  and  because  we  had  established  the 
dominant  transport  and  transformation  processes  for  these  chemicals  in  our 
Phase  II  laboratory  studies.  Since  RDX  was  expected  not  to  transform 
rapidly  in  the  environment,  this  chemical  could  serve  as  a  marker  compound 
to  estimate  the  dilution  process  of  waste  streams  entering  the  Holstoa 
River  from  HAAP. 

A  map  of  the  Ralston.  River  In  the  vicinity  of  the  HAAP  appears  in 
Figure  5.  The  major  input  of  TNT  and  RDX  to  the  Hols  ton  River  comes  from 
the  RDX  lines  1-5  at  Site  A  (river  mile  139.1,  latitude  36*  31’49", 
longitude  82°  38*  17").  We  selected  sices  4.5  km  (Site  B:  river  mile 
136.3,  latitude  36°  31’  30",  longitude  82°  40’  53"),  7.9  km  (Site  C: 
river  mile  134.1,  latitude  36®  28*  52",  longitude  82  41*  02"),  and  33.0  km 
(Site  D:  river  mile  118.5,  latitude  36*  28*  08",  longitude  82*  51*  00") 
downstream  from  Site  A  as  additional  sampling  points.  With  these  selected 
sites  we  expected  the  TNT  concentration  to  show  a  significant  loss  but  the 
RDX  concentration  to  show  only  a  small  loss  by  the  time  wastewaters  dis¬ 
charged  at  Site  A  reached  Site  D.  Therefore,  by  measuring  the  concen¬ 
trations  of  TNT  and  RDX  at  these  four  sites,  we  could  obtain  sufficient 
information  for  the  model  validation. 

We  also  measured  concentrations  of  HMX  (octahydro-1,3, 5,7-tetranitro- 
1,3,5, 7- tetrazocine) ,  SEX  ( 1-ace tyloctahydro-3 , 5 , 7-tr initro-1 ,3,5,7- 

tetrazocine),  and  TAX  ( 1-ace tylhexahydro-3 , 5-d initro-1 ,3, 5-triazine) , 
which  are  discharged  from  HAAP  and  also  appear  at  the  four  collection 
sites.  Although  we  did  not  investigate  the  environmental  fate  of  these 
compounds,  our  analytical  methodology  was  developed  sufficiently  to 
monitor  these  compounds  along  with  TNT  and  RDX.  These  data  may  be  valu¬ 
able  to  the  Army  when  knowledge  about  the  fate  of  these  compounds  is 
desired . 
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B.  Sampling  and  Hydrologic  Bata 

The  collection  of  water  samples  and  hydrologic  data  on  the  Hols  ton 
Elver  was  performed  by  the  U.S.  Geological  Survey  (USGS)  under  the 
direction  of  Sr*  Robert  Llvesey  from  the  Knoxville,  Tennessee  office. 
Since  the  Holston  River  is  approximately  120  meters  wide,  water  samples 
and  hydrologic  data  were  collected  from  boats.  Cables  were  secured  across 
the  river  at  each  sampling  site,  and  markers  were  placed  on  the  cable 
approximately  6  meters  apart.  With  the  boat  attached  to  the  cable  at  the 
marker  position,  the  acquisition  of  samples  and  data  began. 

Water  samples  were  collected  at  each  sampling  location  in  amber 
1-liter  bottles  using  depth-integrating  samplers  (Guy  and  Norman  1970} . 
With  this  procedure,  the  water  was  sampled  as  an  average  of  the  water 
column  at  each  location.  The  bottles  were  labeled  and  stored  in  ice  until 
transferred  to  the  laboratory. 

Water  samples  were  also  collected  at  every  other  sampling  point  for 
turbidity  and  suspended  sediment  analyses.  Turbidity  analyses  were  per¬ 
formed  by  personnel  at  the  University  of  Tennessee  in  Knoxville  and  by 
personnel  at  the  Holston  AAP.  Suspended  sediment  analyses  were  performed 
by  the  USGS  in  Nashville,  Tennessee. 

Sediment  was  collected  from  the  middle  of  the  Holston  River  at  each 
site  along  with  a  water  sample  for  the  determination  of  total  microbial 
counts.  Total  microbial  counts  were  determined  by  Ms.  Kstheryn  Kenyon  of 
the  U.S.  Artsy  Medical  Bioengineering  Research  and  Development  Laboratory 
using  the  serial  dilution  technique  (Bordner  and  Winter  1979).  Sediment 
samples  were  analysed  for  TNT  and  RDX  by  SRI  personnel. 

At  each  sampling  location,  flow  velocity  and  depth  measurements  were 
made  by  USGS  personnel  according  to  standard  procedures  (Buchanan  and 
Somers,  1969).  The  flow  velocity  and  depth  of  the  RDX  lines  1-5  were  also 
determined  by  USGS  personnel.  Flow  velocity  and  depth  determinations  of 
HKX  lines  1-7  and  Arnotc  Creek  were  estimated  by  SRI  personnel. 
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C.  Becenaiaatieit  of  T fff,  8 BX,  KMX,  SEX,  and  TAX 


The  water  staples  collected  from  the  Holscoa  liver  were  analyzed  for 
TNT,  SBX,  HHX  5EX,  and  TAX  using  extraction,  concentration,  and  high- 
performance  liquid  chromatographic  (H1LC)  techniques* 

Approximately  1  liter  of  water  collected  in  an  saber  bottle  with 
teflon  cap  was  filtered  through  a  G,45~is  Millipore  filter  using  vacuum 
filtration.  Samples  of  the  filtrate  (20—25  ml)  were  placed  in  glass  vials 
for  analysis  by  direct  aqueous  injection  into  a  high-performance  liquid 
chromatograph.  The  reraining  filtrate  was  measured  to  the  nearest  5*0  ml 
in  a  graduated  cylinder  and  added  to  a  2.0-liter  separatory  funnel. 

The  eater  sample  was  extracted  with  200  ml  of  dichloromethane 
(Burdick  &  Jackson,  HPLC  grade)  with  vigorous  shaking  for  30  seconds.  The 
extract  was  filtered  through  10-15  gram*  of  anhydrous  sodium  sulfate  using 
vacuum  filtration.  The  extraction  of  the  water  was  repeated  three  times 
with  100-ml  portions  of  dicbloromethane. 

13s*  dried  extracts  were  combined,  transferred  to  a  1-liter  round- 
bottomed  flask,  and  concentrated  to  ~  5  ml  using  a  B-uchi  rotoevapo r« tor 
with  a  water  bath  temperature  of  45°C. 

Tim  extracts  were  transferred  to  8-dram  glass  vials  by  pipet,  the 
flask  rinsed  several  times  with  die hlorome thane ,  and  the  rinses  added  to 
the  vial.  The  vials  were  capped  with  a  Teflon  liners,  wrapped  in  aluminum 
foil,  and  transported  by  air  with  the  water  samples  to  SRI  for  analysis. 

The  extract  was  evaporated  to  dryness  using  a  gentle  stream  of 
nitrogen.  The  residue  was  dissolved  in  0.9  to  1.8  si  of  the  HP1C  mobile 
phase  (65%  water,  35%  methanol),  and  an  internal  standard  (3,5-dinitro- 
toluene)  was  added.  The  resulting  solution  was  analysed  by  EPIC  under  the 
following  conditions: 
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Instrument j 


Spectre-Physics  Model  3500B  Liquid 
Chromatograph 
Column:  Radial  Pak  A~C1S  cartridge  in  a 

Water's  Radial  Compression  Module 

Solvent  program:  Methanol/water  (35/65)  +■  methanol /water 
(65/35)  in  30  ain  linear  gradient 
Plow  rate:  1.2  al/min 

Detection:  UV  <?  254  na 

Retention  Time:  217  a  SEX;  277  s  TAX;  420  s  K(€X;  627  a  RDX; 

1288  s  TKT;  1705  3,5-DNT  (internal  standard) 
Injection  volume:  100-200  pi 

A  typical  chromatographic  profile  appears  in  Figure  6. 

The  concentration  of  each  component  was  determined  by  the  internal 
standard  method  using  peak  area  measurements  determined  by  a  Spectra- 
Physica  "Minigrator"  digital  integrator.  In  some  extracts,  interfering 
peaks  coeluted  with  the  internal  standard.  In  these  cases,  component 
concentrations  were  determined  fey  the  external  standard  method. 

The  concentration  of  each  component  calculated  by  the  peak  area 
determination  was  then  corrected  for  recovery  using  precision  and  accuracy 
data  developed  from  spiked  natural  water  samples. 

Precision  and  accuracy  data  were  obtained  for  all  six  compounds  in  a 
natural  water  (Searsville  Pond,  Stanford,  CA)  over  three  concentration 
levels  including  0.001  mg>_i,  0.025  mgS”1,  and  1.00  tag*""*.  The  percent 
recovery  and  the  standard  deviation  (S.S.D.)  for  each  compound  at  each 
concentration  level  appear  in  Table  1. 

The  percent  recoveries  over  the  three  concentration  levels  were 
averaged  as  shown  in  Table  2  and  these  recovery  factors  were  applied  to 
the  calculation  of  each  munition  component.  With  these  values,  good 
agreement  was  observed  between  those  samples  that  could  he  analysed  by 
direct  aqueous  injection  (no  extraction)  and  the  extract  of  the  same 
sample. 


Tabla  1 

T  Of  K3K1T1CW  WASflWATBJt  CttJHICALS 


rural  wear  froat  Sanravill*  pond,  Stanford 
ilttlw  arandard  deviation. 


u«  Sediment  Analyses 


Sediment  samples  collected  in  glass  jars  were  allowed  to  settle*  and 
the  aqueous  phases  were  decanted  from  the  jars.  An  aliquot  of  the  re¬ 
maining  suspension  was  weighed  and  dried  overnight  at  120*C.  The  before 
and  after  drying  weights  of  sediment  gave  the  water  content  of  the 
sediment.  The  sediment  suspensions  could  therefore  be  handled  without 
altering  the  sediment  properties. 

A  sediment  sample  was  analyzed  by  stirring  approximately  100  g  of 
moist  sediment  with  100  ml  of  ethyl  acetate  In  an  Erletuneyer  flask  for  3 
hours.  The  ethyl  acetate  was  decanted,  an  additional  50  ml  of  solvent  was 
added,  and  the  suspension  was  stirred  for  10  minutes.  The  solvent  was 
decanted,  combined  with  the  first  extract,  and  rotary-evaporated  to 
dryness.  The  residue  was  dissolved  in  a  few  milliliters  of  methanol.  The 
resulting  solution  was  passed  through  a  Clg  Sep-Pak  as  a  cleanup  step  and 
washed  with  a  small  amount  of  methanol,  and  the  eluent  was  evaporated  to 
dryness  with  a  gentle  scream  of  nitrogen.  The  residue  was  dissolved  in  a 
methanol: water  mixture  (50:50)  and  analyzed  by  HPLC  as  described 
previously. 


E.  Confirmatory  Analyses 

The  presence  of  low  levels  of  TNT,  4-&mino~2 ,6-dinitrotoluene  (4-A- 
2,6-DNT),  and  2-amino-4,6~dinitrotoluene  (2-A-4,6-DNT)  found  by  HPLC  was 
confirmed  by  gas  chromatographic  analysis  and  mass  spectronetrlc  analysis 
performed  under  the  following  conditions. 


Instrument : 
Column: 

Temperature : 
Plow  rate: 
Detector: 


Varian  3700  Gas  Chromatograph 

1.8  m  x  2.0  on  glass  column  packed  with  3% 

SP2100  on  80/100  mesh  Chromosorb  W-HP 

180  (4  min  hold)  +  230  @  2*  min”1 

25  ml  min”1  N2 

Electron  capture  (BCD)  or  an  LF»  9000  mass 
spectrometer 
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Table  2 


AVERAGED  PERCENT  RECOVERY  OF  MUNITION  WASTEWATER  CHEMICALS 
FROM  SPIKED  NATURAL  WATER® 


Compound 

SEX 

TAX 

SMX 

RDX 

TNT 


Percent  Recovery  dfc  R.S.D. 
50.4  ±  12.9 
67.3  ±  9.8 
94.2  ±  S.5 

95.9  ±  6.1 

89.9  ±  7.0 


®Natural  water  froa  Searaville  pond ,  Stanford,  CA. 


Gas  chromatography  with  electron  capture  detection  (GC/ECD)  was  based 
on  comparative  analysis  with  the  retention  times  of  authentic  standards. 
Hass  spec trome trie  Identification  was  based  not  only  on  retention  time  but 
also  on  single  ion  monitoring  of  major  Ion  fragments  (m/e  *  210  for  TNT; 
m/e  -  180  for  4-A-2.6-DNT  and  2-A-4,6-DNT) . 

Samples  were  collected  from  the  HPLC  according  to  their  elution 
volume ,  concentrated  to  dryness  using  nitrogen  gas,  and  dissolved  In  a 
small  amount  of  methanol  before  gas  chromatographic  analysis. 

F.  Sate  and  Equilibrium  Constants 

The  photochemical  rate  constant  for  TNT  was  measured  In  North  and 
South  Fork  Holston  Elver  water  obtained  upstream  from  the  HAAP  discharge 
lines.  TNT  was  spiked  into  the  water,  and  photolysis  experiments  in 
sunlight  were  performed  under  several  conditions—*  in  filtered  sterilised 
water  in  borosillcace  tubes  and  in  10  cm  x  200  cm  open  crystallizing 
dishes  and  in  unfiltered  water  containing  suspended  sediment  in  open 
dishes.  The  TNT  concentrations  were  measured  as  a  function  of  tine,  and 
the  rate  constant  was  calculated  from  the  integrated  form  of  Equation  1 
shown  in  Equation  8, 

m  -  kpt  (8) 


where  CQ  and  Ct  are  the  TNT  concentrations  at  times  zero  and  t. 
Regression  of  In  C  versus  t  gives  the  slope  kp. 

Pseudo-first-order  rate  constants  for  the  photolysis  of  RDX  and  the 
blotransformatioa  of  TNT,  as  well  as  the  TNT  and  KDX  sorption  partition 
coefficients,  were  obtained  from  our  Phase  II  report  (Spaaggord  et  al« , 
1980b).  A  second-order  rate  constant  for  the  biotransformation  of  TNT  was 
obtained  using  Equation  5. 
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IV  RESULTS 


The  field  sampling  studies  of  the  Holston  River  took  place  on  April 
22  and  23,  1981.  Weather  conditions  were  hazy  and  sunny  on  April  22  and 
rainy  oa  April  23.  Rains  from  earlier  in  the  week  in  West  Virginia 
created  above  average  flows  (runoff)  in  the  North  Fork  of  the  Holston 
River;  however,  the  river  receded  during  the  collection  of  data. 

On  April  22,  the  total  flow  contributions  from  the  North  and  South 
Forks  of  the  Holston  River  were  71.5  m3s”1  and  22.2  m^s"1  respectively. 
On  April  23,  flows  were  73.5  m3s_1  and  21.6  ta3s_1,  respectively.  These 

values  were  obtained  from  upstream  gaging  stations  on  each  fork.  Flow 

measurements  downstream  at  Surgoinsville  gaging  station  (Site  D)  were 
120.8  m3  s'*1  on  April  22  and  89.9  a3  s'1  on  April  23. 

A.  Hydrological  Data 

1.  Streambed  Contour  and  River  Flow  Velocity 

From  the  depth  measurements  at  each  sampling  location  (Table  3),  we 
constructed  a  streambed  contour  for  each  of  the  four  sites  investigated. 
Contours  for  sites  A,  B,  C,  and  D  appear  in  Figure  7.  At  each  section, 
the  cross-sectional  area  was  calculated  (distance  between  sections  times 
depth)  and  summed  across  the  width  of  the  river  to  yield  a  total  cross- 

sectional  area.  The  total  cross-sectional  area  was  divided  by  the  total 

summed  flow  rate  to  yield  the  river  flow  velocity.  These  data  appear  in 
Table  4.  The  flow  velocity  was  used  to  estimate  the  residence  time  of  the 
water  in  each  compartment. 

The  flow  rates  at  sites  A,  B,  and  C  vary  within  42,  which  may  be  due 
to  the  natural  transportation  of  water  into  and  out  of  each  site  or  the 
accuracy  of  the  measurement.  However,  a  large  reduction  in  total  flow  was 
noted  between  Sites  C  and  D  (17X).  Seepage  may  be  one  factor  contributing 
to  the  loss  of  water;  however,  no  causes  for  this  loss  were  Identified. 

2.  Water  Quality  Parameters 

Water  quality  parameters  including  temperature,  pH,  dissolved  oxygen, 
and  conductivity  were  measured  at  each  sampling  station.  These  data 

appear  in  Table  5.  These  data  indicate  that  the  Holston  River  ia  in  good 
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Table  3 


DEPTH  AND  PLOW  MEASUREMENTS  ACROSS  THE  HOUSTON  RIVER 


AT  POOR 

SAMPLING 

LOCATIONS 

Section 

Station* 

Time 

Depth 

(feet) 

ps* 

(ft3rx 

A-  1 

25 

9:20 

5.7 

63.9 

A-  2 

40 

9:30 

6.0 

132 

A-  3 

55 

9:35 

6.0 

139 

A-  4 

70 

9:45 

5.3 

161 

A-  5 

90 

9:50 

6.1 

238 

A-  6 

110 

9:55 

5.9 

203 

A-  7 

130 

10:00 

6.0 

172 

A-  S 

150 

10:10 

5.7 

182 

A-  9 

170 

10;15 

5.5 

167 

A-10 

190 

10:20 

5.7 

173 

A-ll 

210 

10:27 

5.5 

190 

A-12 

230 

10:32 

5.7 

189 

A-13 

250 

10:45 

6.3 

214 

A-14 

270 

10:52 

5.7 

201 

A-15 

290 

11:00 

5.2 

180 

A-16 

310 

11:05 

5.1 

199 

A“i7 

335 

11:09 

5.6 

162 

A-18 

360 

11:18 

5.3 

169 

A-19 

385 

11:25 

5.0 

140 

A-20 

410 

11:30 

4.5 

190 

A-21 

450 

11:35 

4.0 

69.6 

Total  Flow 

3524.5 

aec) 


Table  3  (Continued) 


Table  3  (Continued) 


Be 


Site  C  Section  Station 


Time 


LA-7934-6 


FIGURE  7  STREAMBSD  CONTOURS  OF  EACH  SAMPLING  LOCATION 
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ecological  condition.  Some  divergence  of  the  data — auch  as  that  on  con¬ 
ductivity  at  Sites  A,  8,  and  C— •indicates  that  the  river  ia  not  well  aimed 
within  the  boundaries  of  HAAP. 

3 .  Suspended  Sediments  and  Turbidity 

The  suspended  sediment  and  turbidity  were  measured  at  several 
stations  at  each  site.  These  data  appear  in  Table  6.  The  suspended 
sediments  were  found  to  decreese  from  approximately  126  mg!-”1  at  Site  k  to 
22  at  Site  0.  It  is  expected  that  the  bottom  sediment  layer  will 
increase  under  these  conditions.  Because  both  TOT  and  RDX  sorbed  to 
sediment  weakly,  we  expect  only  small  losses  of  these  chemicals  from  the 
water  to  the  sediment.  Turbidity  also  decreased  from  Site  &  to  Site  D  in 
a  fashion  similar  to  the  sedimentation. 


The  microbial  population  was  determined  from  one  section  at  each 
sampling  site.  These  results  appear  in  Table  7.  The  value*  (10s  to  10** 
cells  si"1)  are  low  compared  to  the  value  of  106  cells  ml"1  usually  found 
in  a  running  stream.  Higher  values  ware  expected  because  of  waste  from 
upstream  industries  and  the  high  sediment  loadings.  However,  the  observed 
values  were  used  in  the  model,  end  therefore  biotrensformatlon  was 
expected  to  contribute  little  to  the  overall  loss  of  t«T  or  RDX. 


B.  Input  Concentrations  of  Chemicals  from  BMP  Westeetreama 

Thrae  was tea cream*  provide  chamlcal  input  into  tha  Holeton  River  from 
HASP.  These  westestreams  are  the  HMX  lines  1-7,  RDX  lines  1-5,  and  Amott 
Creek.  RDX  lines  1-5  were  believed  to  be  the  major  Input  of  TOT  and 
RDX.  The  HMX  lines  1-7  enter  the  Holston  River  approximately  600  m  above 
the  RDX  lines,  while  Amott  Creek  enters  the  Holston  River  approximately  4 
fas  downstream  from  the  RDX  lints.  A  map  daaignatiog  the  input  streams  Is 
shown  in  figure  8. 

The  RDX  input  stream  travels  through  an  underground  culvert  and 
surfaces  approximately  13  m  from  the  Holston  River,  where  it  mixes  with 
cooling  water.  The  flow  of  the  scream  prior  to  miring  with  the  cooling 
water  was  estimated  et  0.07  m3  e"1  (2.63  cfs). 
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Table  6 

SUSPENDED  SEDIMENT  AND  TURBIDITY  MEASUREMENTS 
AT  SITES  A,  B,  C,  AND  D 


S.S.a$j£g, 


2 

B 

KB 

8 

10 

12 

14 

16 

18 

Suspended  Sediaents  (tag  l~l) 

121 

126 

126 

124 

116 

109 

103 

75 

88 

236 

68 

99 

98 

93 

90 

92 

64 

81 

67 

91 

90 

94 

90 

87 

37 

79 

72 

22 

29 

31 

32 

35 

38 

34 

44 

43 

Turbidity  (tag  t~  ) 


14 

12 

18 

15 

12 

11 

12 

8 

6 

3 

21 

17 

18 

23 

26 

21 

21 

22 

21 

23 

10 

11 

12 

11 

12 

12 

11 

8 

8 

2 
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Table  7 


MICROBIAL  POPULATIONS  MEASURED  IN  THE  HQLSTON  RIVER 
AT  EACH  SAMPLING  SITE 


Sample 

Dilution 

Raw  Counts a 

Avg  ral”^ 

CPU  ml“lc 

A-7 

10° 

228,  253,  265 

2487 

10' 1 

63,  98,  39 

667 

— 

10~* 

5,  6,  2 

43 

4.5  x  10* 

io~s 

0,  0,  0 

— 

— 

10~“ 

~0 

o 

o 

o 

B-13 

10 

TNTC  x  3 

— 

— 

10“ 1 

104,  90,  102 

987 

— 

10“  * 

9,  16,  11 

120 

1.1  X  10* 

10”  3 

1,  1,  1 

10 

MHH 

10““ 

0,  0,  0 

— 

— 

C-ll 

10° 

441,  456,  351 

4160 

— 

10-' 

56,  60,  48 

547 

— 

10“a 

3,  11,  7 

70 

4.2  x  10* 

10”  3 

o 

H* 

O 

3 

— 

10”“ 

0,  0,  0 

— 

— 

D-10 

10° 

134,  100,  104 

1127 

10' 1 

11,  4,  11 

87 

— 

10”* 

1.  1,  2 

13 

1.6  x  10* 

10”* 

0,  1,  0 

3 

— 

10~“ 

0,  0,  0 

— 

— 

& 

Colonies  per  plate. 

k Average  number  of  colonies  per  milliter  of  sample. 
cColcny  forming  units  per  milliter  of  Holston  River  Water. 
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The  HHX  lines  1-7  are  also  underground  until  approximately  450  a  from 
the  Holaton  Elver,  where  che  rater  surfaces,  mixes  with  cooling  water,  and 
travels  toward  the  river.  The  flow  of  the  HKX  lines  1-7  was  not  measured 
but  estimated  at  0.07  a3  s'*1,  then  mixed  with  cooling  water,  a  stream 
approximately  1.5  ra  x  0.3  m  is  produced  with  an  estimated  velocity  of  0.45 
m  a**1. 

Arnott  Creek  la  approximately  7.6  a  by  0*45  m  and  travels  approxi¬ 
mately  900  a  to  the  Holaton  Elver.  Wastes  emptied  into  Arnott  Creek  are 
usually  from  building  washout  operations.  The  flow  rate  of  the  creek  was 
estimated  at  0.83  m3  a""1. 

The  concentrations  of  TNT  and  RDX  were  determined  for  each 
wastestream  by  HPLC  from  samples  taken  midway  between  the  plant  and  the 
Holaton  River  in  Arnott  Creek  and  from  the  RDX  and  HMX  lines  prior  to 
sunlight  exposure.  The  input  concentrations  and  mass  calculations  appear 
in  Table  8. 

Table  8 

TNT  AND  RDX  DISCHARGES  INTO  THE  H0LST0N  RIVER 


Wastestream 

Discharge 
Rate 
(a3  s”1) 

_  TNT 

mg  l”1  «g  8-1 

-  RDX 

mg  l”l  mg  s'* 

HMX  lines  1-7 

0.07 

0.23 

16.1 

8.37 

586 

BDX  lines  1-5 

0.07 

0.64 

44.8 

2.49 

181 

Arnott  Creek 

0.83 

0.002 

1.7 

0.18 

149 

C.  Rate  and  Equilibrium  Constanta 

1.  TNT  Photolysis  Rate  Constant 

To  determine  a  photolysis  rate  constant  for  the  TNT  model  simulation, 
we  measured  the  sunlight  photolysis  rate  constants  of  TNT  in  natural 
waters  in  several  types  of  reaction  cells  and  in  the  presence  of  suspended 
sediments  to  determine  how  these  factors  may  affect  photolysis  rates. 
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To  simulate  an  environmental  situation  In  which  sunlight  irradiates  a 
quiescent  water  body,  ve  chose  s  flat-bottomed  crystallising  dish  5  cm 
deep  and  10  cm  In  diameter.  Table  9  lists  data  on  the  photolysis  of  TNT 
in  waters  from  the  South  Fork  of  the  Holston  River  in  open  dishes  and  in  a 
pyrex  tube.  The  two  dish  experiments  used  filtered  South  Fork  water  in 
one  dish  and  unfiltered  water  in  the  other  dish,  in  which  the  sediment 
(collected  with  the  sample)  was  kept  in  suspension  by  stirring.  The 
photolyses  were  conducted  on  a  clear  sunny  July  day  between  10:30  AH  and 
1:30  FM  (high  sun  angle);  therefore,  the  solutions  In  the  dishes  were 
irradiated  through  the  surface  of  the  water  without  any  filtering  of  the 
light  through  glass. 


Table  9 


PHOTOLYSIS  OF  0.61  ppm  TNT  IN  NATURAL  WATERS  IN  PRESENCE 
OF  PARTICULATES  AND  IN  SEVERAL  TYPES  OF  REACTION  CELLS® 


Solutions  and  Conditions 

Cell 

Measured  First-Order 

Rate  Constants,  kp  day”1 

Holston  River  water 

(South  Fork),  0.2-) 

Tube 

89 

filtered 

Dish 

26 

Holston  River  water 

(South  Fork), 

unfiltered  and  sedi¬ 
ment  suspended  by  stirring 

Dish 

18 

®Photolysed  in  July  between  10:30  A.M.-  1:30  P.M. 
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The  data  show  that  photolysis  of  TNT  in  filtered  natural  waters 
occurs  about  3.4  times  faster  in  the  tube  than  in  the  dish.  The 

difference  In  these  photolysis  rate  constants  is  mainly  because  of  the 
greater  light  flux  passing  through  the  tube  (360*  exposure  to  light)  than 
through  the  dish  (light  enters  only  from  open  top  of  dish). 

The  data  in  Table  9  for  the  Holston  Elver  water  alao  indicate  that  in 
unfiltered  waters  containing  sediment,  the  photolysis  of  TNT  will  be  about 
50Z  slower  chan  that  in  filtered  natural  water.  The  slower  photolysis 
rate  is  due  both  to  light  scattering  and  to  light  absorption  by  the 

particulates.  Our  understanding  of  the  effects  of  particulates  on 
photolysis  rates  in  natural  water  Is  incomplete,  but  obviously  such 

effects  will  be  impoitant  In  deep  or  turbid  systems.  Therefore,  the 

maximum  photolysis  race  for  TNT  in  the  Holston  Elver  is  expected  to  be  18 
day”1  in  July.  Using  the  data  of  Zepp  and  Cline  (1977)  for  light 
intensity  in  late  spring  (to  correspond  to  field  sampling  dates),  the 

photolysis  rate  constant  reduces  to  16  day'*1. 

In  our  laboratory  and  sunlight  photolysis  studies  of  TNT,  the 

experiments  have  been  performed  under  conditions  in  which  the  solutions 
are  optically  "thin";  such  solutions  weakly  absorb  light  with  most 

incident  light  passing  through.  In  deep  aquatic  system  of  some  depth, 
however,  light  available  for  initiating  reactions  will  decrease  with 
increasing  depth.  Therefore,  the  photolysis  rate  of  TNT  will  decrease 
with  depth,  and  obviously  no  photolysis  will  occur  below  the  level  where 
all  light  is  absorbed.  Solutions  In  which  most  of  the  incident  light  is 
absorbed  are  referred  to  as  optically  "thick”;  it  should  be  noted  that 
optically  "thick”  solutions  may  result  from  a  long  pathleagth  or  a  high 
absorption  coefficient  of  the  solution  or  both. 

As  an  example  of  light  penetration  into  a  natural  water,  we  have 

calculated  the  depth  of  Holston  Elver  water  at  which  99X  of  the  Incident 
light  at  the  surface  is  absorbed.  This  situation  occurs  when  the 
absorptivity  of  the  water  exceeds  2.0:  that  is  when 


aK  x  0  >2 


(9) 
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where  Is  the  absorptivity  at  wavelength  X  and  B  is  the  depth  In  ceati- 
meter a.  Table  10  lists  the  depths  for  several  wavelengths  in  the  solar 
spectrum* 


Table  10 

DEPTHS  AT  WHICH  99X  OP  LIGHT  IS  ABSORBED 
BY  HOLSTOH  RIVER  WATER* 


Wavelength 

X^(nm) 

Absorbance 
a^(cm~* ) 

Depth 

(cm) 

300 

0.045 

45 

350 

0.021 

96 

400 

0.011 

180 

450 

0.006 

330 

500 

0.004 

530 

As  these  data  show,  the  longer  (and  less  energetic)  wavelengths 
penetrate  deeper  into  the  Hoiston  River  water  than  do  the  shorter  and  more 
energetic  wavelengths. 

If  sunlight  were  simply  absorbed  by  the  natural  water  without 
promoting  the  photolysis  of  TNT,  the  effect  on  the  direct  photolysis  rate 
constant  of  TNT  would  be  a  straightforward  light  screening  phenomenon  that 
reduces  the  photolysis  rate  of  TNT  (Mill  et  al.,  1981;  Zepp  and  Cline, 
1977).  However,  since  light  absorption  by  the  natural  water  is  respon¬ 
sible  for  the  transformation  of  TNT,  the  estimation  of  the  TNT  loss  rate 
due  to  photolytlc  processes  Is  more  complex. 

As  defined  and  developed  in  Part  II-C  of  this  report  ,  Photo¬ 
chemical  Studies  of  TNT  In  Natural  Waters,  the  reaction  quantum  yield 
for  indirect  photolysis  of  TNT  is  given  by  Equation  10.  Rearranging  this 
equation  gives  the  rate  of  loss  of  TNT: 


*Sample  taken  in  April,  1981 
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“*  fit  JL 

-  k  [Tin  ]  -  4.1  [1-10  X  ]  {TNT]  (10) 

tl  t  p 

As  shown  by  the  data  In  Table  11,  the  calculated  0  values  at  313  na  and 
366  as  were  roughly  the  same  at  the  same  Initial  TNT  concentrations. 
Unfortunately,  we  could  not  perform  enough  experiments  to  determine  at 
what  longer  wavelength  the  6  value  would  decrease  such  that  photolysis 
rate  of  TKT  would  reach  a  limiting  value  (that  Is,  at  some  less  energetic 
wavelength,  photolysis  of  TNT  would  be  less  efficient  or  not  occur  due  to 
insufficient  energy).  More  studies  need  to  be  performed  before  this 
approach  tc  estimating  TNT  photolysis  loss  rates  can  be  verified. 


Table  11 

QUANTUM  YIELDS  FOR  TNT  MEASURED  IN  SEARS VILLE 
POND  WATER  AT  313  NM  AND  366  NM 

fTST]o  l(rna)  4» 

2.79  x  lO"5  31 3a  8.8  x  10“ 

2.79  x  10~«  366b  5.0  x  10”** 

An  alternative  procedure  for  calculating  the  photolysis  rate  constant 
of  TNT  in  Holston  River  water  (HRW;  HRW*  indicates  an  excited  state) 
considers  the  following  model  scheme: 

HRW  -»  HRW* 

HRW*  +  TNT  TNT  products  +  HRW 

HEW*  +  HEW  kd 


aI.  -  1.37  x  IQ6  Einsteins  Jt1  sec1  at  313  ns 
“  9*85  *  106  Einsteins  A1  sec1  at  366  no 
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where  ka  is  the  rate  constant  for  light  absorption  by  the  natural  water , 
kr  is  the  rate  constant  for  production  of  products  from  reaction  of  TUT 
and  HEW*,  and  kg  is  the  rate  constant  for  decay  of  KHW*.  The  total  rate 
of  loss  of  TUT  is  defined  by  Equation  11, 


— d(YST) 


dt 


k  k 

-j—  (TNT) 
Kd 


kpE(Ttfr> 


<u) 


where  kpE  is  a  first-order  rate  constant  for  loss  of  TNT. 


Laboratory  experiments  that  use  1-cm  diameter  tubes  for  photolysis 

give  values  of  fcpE  for  optically  thin  natural  water  solutions  (k*g)  .  The 

rate  constant  of  interest  for  modeling  the  fate  of  T8T  in  deep  aquatic 

systems  is,  however,  that  for  the  optically  thick  solution  (  k7?  )•  The 

pE 

relationship  between  the  two  rate  constants  for  midday  sunlit  conditions 
cen  be  shown  to  be  dependent  only  on  the  depth,  D,  and  the  values  of 
the  water  body  according  to  Equation  12, 


,TN 

_El_ 

pE 


•  2.3  D 


(12) 


where  ax  is  the  midday  sunlight  intensity  corrected  for  pathlength  in  a 
water  body  and  is  the  midday  sunlight  intensity  incident  on  the  water¬ 
body.  Values  of  Z k  and  W^,  listed  by  Zepp  and  Cline  (1977)  for  spring 
season  dates  at  40°  latitude  and  290  am  to  500  nm,  were  used  at  tr^  values 
measured  for  the  Holston  River  water  (South  Fork),  to  calculate  the  value 
of  Equation  13. 


kTN 

»  0.0220  (13) 
kps 

This  equation  indicates  that  for  the  spring  season  at  midday  and  40* 

latitude  the  two  first-order  rate  constants  are  equal  at  a  depth  of  45  cm; 

TK  TH 

in  deeper  waters  k  _  will  be  smaller  than  k  _  . 
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fox  TOT  for  use  in  a  face  model  of  the 


Tl? 

To  estimate  a  value  of  k  „ 

pS 

Ho  Is  ton  River,  we  assumed  an  average  river  depth  of  180  cm  <6  ft),  used  a 
—1  TH 

value  of  16  day  for  k  _  as  determined  previously,  and  solved  Equation 

13  for  k  _  to  yield  an  environmental  photolysis  rata  constant  of  4.0 
.  Ps 

day  at  midday. 

2.  Other  Race  and  Equilibrium  Co as tents 

llie  remaining  rate  and  equilibrium  constants  along  with  the  1ST  rate 
constant  used  in  this  study  appear  in  Table  12.  The  sediment  sorption 
partition  coefficients  were  obtained  from  the  Phase  II  study  (Spanggord  et 
al.,  1980b).  The  blotransfaraation  rate  constant  was  also  taken  from  the 
Phase  II  study;  however,  due  to  the  observed  low  microbial  counts,  this 
constant  was  small  compared  to  that  for  photolysis.  The  photolysis  rate 
constant  for  RDX  was  also  obtained  from  the  Phase  XX  study. 


Table  12 


RATE  AND  EQUILIBRIUM  CONSTANTS  FOR  TNT  AND  RDX 
USED  IN  THE  SRI  MODEL 


Chemical 


Partition  Photolysis 
Coefficient  (day-1) 


Bio transformation 
(day-1) 


RDX 

TOT 


5.0 

0.1 

— 

50 

16.4 

0.03 

D.  Simulation  of  the  One-Dimensional  Compartmentalized  River  System 

For  the  one— dimensional  compartmentalised  river  system,  the  Holston 
River  was  divided  into  67  compartments  of  equal  length  over  the  33  km  of 
the  studied  river  reach.  The  reach  was  characterised  as  four  different 
segments  with  each  segment  described  by  the  geometric  data  surveyed  at  the 
sampling  site*  These  characterizations  are  shown  in  Table  13. 


table  13 


GEOMETRICAL  IPW  DATA  FOR  TEE  ONE-DIMENSIONAL 
CGMPAETKENTALI ZW  Rim  SYSTEM 


Reach  Compartments 


Volume  Settling  Rate  Scouring  Rate 
(JO  (%  day*1)  (X  day*1) 


1  1  to  9 

2  10  to  18 

3  19  to  40 

4  41  to  67 


1.0  x  108  0.82 

1.3  x  IQS  0.64 

0.8  x  10®  0.64 

1.4  x  10s  0.S0 


2.0  x  10“** 
2.0  x  10~* 
2.0  x  10“* 
1.4  x  1Q~* 


The  settling  rate  in  each  segment  was  estimated  from  the  suspended 
sediment  measurements  made  at  each  site.  We  assumed  that  the  sediments  were 
in  equilibrium  during  the  simulation  period.  The  scouring  rate  was  estimated 
by  assuming  that  IX  of  the  suspended  solids  were  returned  to  the  aqueous  phase 
from  the  bottom  sediment. 


RDX  was  selected  for  the  initial  investigation  of  the  one-dimensional 
system  because  RDX  is  a  relatively  conservative  chemical,  and  therefore,  the 
mixing  and  dilution  properties  of  the  Eolston  River  could  be  evaluated.  The 
RDX  loadings  from  the  HHX  lines,  RDX  lines,  and  Arnott  Creek  were  entered  into 
compartments  1,  2,  and  7,  respectively.  The  concentration-time  (distance) 
dependent  profile  of  RDX  predicted  by  the  one-dimensional  system  appears  in 
Figure  9. 

At  Site  A,  the  predicted  RDX  concentration  is  8.0  ppb.  Ths  actual  values 
found  ranged  from  212  ppb  (North  Bank)  to  less  than  1.0  ppb  on  the  South  Bank 
(see  Section  1V-G) .  The  actual  values  found  across  the  width  of  the  river 
indicate  that  the  Rolston  River  Is  not  well  mixed.  Since  the  one-dimensional 
system  assumes  a  well-nixed  reaction  vessel,  the  application  of  the  model 
before  the  river  becomes  homogeneous  will  lead  to  erroneous  prediction*.  The 
predicted  and  observed  values  of  RDX  at  each  sampling  site  appear  in  Table  14. 
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COMPARTMENT  NUMBER 


1 _ L. . -J . . L 

0  5  10  IS 

DISTANCE 


20 

km 


25  SO 

LA-7S34-8 


FIGURE  9  RDX  CONCC NT ft  AT! ON  PROFILE  SN  THE  HOLSTON  RIVER 
PREDICTS  0  BY  ONE-OI SSENSIONAL  MODEL 


WVZ’trrbgB  n 


Table  14 


ONE-DIMENSIONAL  COMPARTMENTALIZED  RIVER  SYSTEM 
PREDICTIONS  AND  FOUND  CONCENTRATIONS  OF 
RDX  IN  TEE  HOLSTON  RIVER 


Site 

Predicted  Concentration ppb) 

Found  Concentration  ppb) 

A 

8.0 

North  Bank 

212. 

Middle 

1.27 

South  Bank 

0.5 

B 

9.6 

7.3 

12.3 

4.1 

C 

9.5 

12.5 

10.3 

4.6 

D 

9.3 

j.7 

7.1 

7.6 

From  che  observed  concentrations  of  RDX,  It  is  clear  that  homogeneity 
in  the  Holston  River  Is  eventually  achieved  somewhere  between  Sites  C  and 
D.  At  Site  U»  where  homogeneity  is  achieved,  the  predicted  value  is  in 
good  agreement  with  the  observed  value. 

B.  Simulation  of  the  Two-Dimens tonal  Compartmentalized  River  System 

Since  the  one-dimens ional  compartmentalized  river  system  could  not 
adequately  describe  the  concentration  profile  of  RDX  near  the  discharge 
sites,  we  employed  a  two-dimensional  compartmentalized  river  system  that 
considers  both  lateral  and  directional  flows  to  describe  the  loss  and 
movement  of  TNT  and  RDX  in  the  Holston  River. 

1.  Comparta&ntalization  and  Lateral  Flow 

The  Holston  River  was  compartmentalized  as  shown  in  Figure  10,  where 
the  left-hand  compartments  represent  the  North  ftanV  and  the  right-hand 
compartments  repres  at  the  South  Bank.  In  the  tw**  dimensional  compart¬ 
mentalized  river  system,  the  chemicals  Introduced  .  to  compartment  1  are 
allowed  to  move  laterally  into  compartment  2  as  well  as  in  the  direction 
of  the  river  flow  (compartment  4).  The  HAAF  discharge  lines  enter  the 
two-dimensional  compartmentalized  river  system  at  compartments  1  (HMX 
lines  1-7),  4  (RDX  lines  1-5),  and  16  (Arnott  Creek).  The  lateral  flow 
was  estimated  from  the  dye  studies  of  Sullivan  et  al.  (1977)  who 
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determined  that  discharges  from  the  RDX  lines  1-5  would  not  become 
homogeneous  until  10  km  from  the  discharge  point  (with  an  irta cream  flow  of 
223  a3  s'*1).  Assuming  an  average  flow  velocity  of  0.43  m  s”1 ,  we 

calculated  it  would  take  6  hours  for  water  to  traverse  10  km  and  therefore 
6  hours  to  traverse  the  width  of  the  river.  Therefore,  for  Sullivan’s 
study,  we  can  estimate  the  lateral  flows  at  Sites  A,  B,  and  C  to  he  4.5  x 
i0~3  m3  s'*1.  During  the  field  study,  the  Instreea  flow  was  only  93 
at3  g"*!,  or  approximately  half  of  that  observed  by  Sullivan.  We  therefore 
assumed  that  our  lateral  flow  velocities  would  be  half  of  those  observed 
by  Sullivan  and  that  homogeneity  of  the  Holston  River  would  not  be 

achieved  until  20  ka  downstream  from  the  discharge  point.  After  20  km, 
the  one-dimensional  compartmentalized  river  system  wav  employed  for  the 
remainder  of  the  river  reach. 

2*  Cross-Sectional  Area  of  Each  Compartment 

To  determine  the  cross-sectional  area  of  each  compartment,  we  plotted 
the  cumulative  flow  against  the  cumulative  cross-sectional  area  from  the 
data  collected  by  the  0SCS  at  each  site  (Table  1).  This  plot  appears  in 
Figure  11.  We  arbitrarily  divided  the  river  into  three  lateral 

compartments  with  flows  of  16.7,  25.0,  and  53  m3  a”1  in  the  transect 

direction.  We  then  computed  the  cross-sectional  area  and  flow  velocity 
for  each  compartment;  these  data  are  listed  in  Table  15.  These  data  show 
a  greater  flow  in  the  center  of  the  river,  Which  is  consistent  with  the 
observed  results. 


Table  15 


CROSS-SECTIONAL  AREAS  AND  FLOW  VELOCITIES  OF 
LATERAL  COMPARTMENTS  AT  SITES  A,  B*  AND  C 


Site 

Flow 

iw  cb. 

X-Area 

<m2) 

A 

Velocity 
(«  •-*> 

X-Area 

-  j  — -  w 

B 

Velocity 
<m  s”1) 

X-Area 

C 

Velocity 

(a  *~1) 

16.7 

36.7 

0.46 

55 

0.30 

27.5 

0.61 

25.0 

50.5 

0.50 

60 

0.42 

41.3 

0.61 

53.0 

119.4 

0.44 

147 

0.36 

99.2 

0.53 

3.  RDX  Simulation 

The  results  for  the  simulation  of  RDX  concentrations  are  shown 
graphically  for  the  North  Bank,  Center,  and  South  Bank  of  the  Rolston 
River  in  Figure  12.  The  profiles  indicate  that  the  North  Bank  undergoes 
continuous  dilution  until  homogeneity  in  the  river  Is  achieved  at  20  km. 
The  concentration  of  RDX  In  the  center  of  the  river  increases ,  exceeds  a 
downstream  homogeneous  concentration,  then  tapers  off  to  a  constant  value 
(homogeneity).  The  South  Bank  RDX  concentration  Increases  continually 
until  homogeneity  ie  achieved.  The  transport  (sorption)  and  trans¬ 
formation  (photolysis)  processes  have  little  effect  on  reducing  the  RDX 
concentration  over  the  studied  river  reach.  For  example,  the  RDX  concen¬ 
tration  at  20-k»,  9.2  ppb,  will  decrease  to  9.0  ppb  at  Site  D.  We  expect 
only  5%  of  the  RDX  to  be  transformed  over  33  km.  Therefore,  In  the  case 
of  RDX,  the  two-dimensional  compartmentalised  river  system  describes  the 
dilution  and  mixing  properties  of  the  Holston  River. 

F.  TNT  Simulation  Employing  a  Two-Dimensional  Compartmentalised 

liver  "System .  *  '  ""rf . r 

Using  the  Input  concentration  of  Tiff  found  at  Site  A-l,  we  ran  the 
TNT  simulation  two-dimensional  compartmentalised  river  system  using 
photochemical  rate  constants  of  A  and  16  day"1.  The  concentratlon-tlme- 
of-day-dependent  curves  for  each  rate  constant  appear  in  Figure  13  and  14 
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foe  the  North  Bank  of  Sice  A  and  the  North,  Center,  and  South  Backs  of 
Sites  B  and  C.  These  curves  readily  show  the  TNT  concentration  dependence 
on  the  photochemical  rate  constant.  Although  the  value  of  4  day-1  was 
used  in  this  study  {as  determined  from  the  laboratory  studies),  this  value 
was  generated  from  highly  turbid  waters  due  to  rains  and  runoff  In  the 
area.  We  expect,  therefore,  that  under  less  turbid  conditions,  the 
photochemical  rate  constant  may  rise  to  16  day"1  or  above  in  the  Holston 
River. 

Using  the  photochemical  rate  constant  of  4  day"*1 ,  the  predicted 
profile  of  TNT  concentration  as  a  function  of  distance  from  the  discharge 
point  is  shown  in  Figure  15  for  the  North,  Center,  and  South  portions  of 
the  Holston  Elver  at  noontime.  The  curves  are  highly  divergent  at  Site  A 
due  to  poor  mixing,  and  follow  a  similar  pattern  to  that  of  KDX  (Figure 
12).  The  curves  converge  at  20  km  (homogeneity),  then  decrease  as  a 
function  of  the  photochemical  rate  constant.  We  expect  that  TNT 
concentrations  near  the  plant  boundary  (8  km)  are  primarily  controlled  by 
dilution  under  the  conditions  of  this  study. 

G.  Monitoring  Data  From  The  Bolaton  River 

The  concentrations  of  munitions-related  pollutants  being  discharged 
into  the  Holston  Elver  were  determined  by  HPLC  as  described  in  Section 
1X1.  The  components  were  identified  by  their  capacity  factors  as  compared 
to  authentic  standards.  The  chromatographic  profiles  for  all  the  analysed 
samples  appear  in  Appendix  1.  We  observed  that  the  chromatographic 
profiles  became  complex  close  to  the  South  Bank.  This  probably  results 
from  wastes  discharged  from  upstream  industries  located  on  the  South  Fork 
of  the  Holston  Elver.  In  some  chroaatograes  the  waste  products  coeluted 
with  our  internal  standard.  We  therefore  calculated  the  concentration  of 
the  munitions-related  pollutants  by  both  internal  and  external  standard 
methods.  The  observed  concentrations  of  SEX,  TAX,  HHX,  RBX,  and  TNT  at 
each  site  and  section  appear  in  Table  16. 

The  data  indicate  that  at  Site  A,  moat  of  the  discharge  is  hugging 
the  North  Bank  with  little  dispersion  into  the  river.  At  Sites  B  and  C, 
the  lateral  transport  of  chemical  becomes  evident,  and  at  Site  D,  complete 
mixing  of  the  river  has  been  achieved.  The  dispersion  phenomenon  is 
depicted  in  Figure  16  for  EDX. 
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FIGURE  15  TNT  CONCENTRATION  AT  NOON  AS  A  FUNCTION  OF  DISTANCE 
FROM  THE  DISCHARGE  POINT  FOR  THE  NORTH,  CENTRAL, 

AND  SOUTH  PORTIONS  OF  THE  HOUSTON  RIVER 
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Hie  data  also  show  thee  little  Wf  and  SIX  were  being  discharged 
daring  the  day  of  sample  collection  (April  22).  At  all  four  sites,  many 
of  the  TUT  analysis  values  were  at  background  levels*  This  background 
concentration  averaged  to  0.83  ppb.  Random  GC-SCO  analyses  of  sample 
extracts  confirmed  the  presence  of  TUT,  hot  these  results  could  not  be 
confirmed  by  mass  spectrometry. 

Sue  to  the  rapid  transformation  of  THT  in  the  environment  mod  the 
large  dilation  factors  caused  by  heavy  flow,  THT  could  he  detected  only  at 
a  few  locations.  The  paucity  of  analytical  data  was  sufficient  to  hinder 
rigorous  validation  of  the  SRI  Model  using  THT. 

H.  Sediment  Analyses 

The  sediment  extracts  were  analyzed  by  HPLC  as  previously  des¬ 
cribed.  Very  large  peaks  eluted  with  retention  volumes  similiar  to  those 
of  SEX,  TAX,  and  HMX,  and  the  presence  or  absence  of  these  compounds  was 
not  confirmed.  RDX  was  identified  in  the  sediments  along  with  4-eaino- 
2, 4-dinitro toluene  and  2 -amino-4 , 6—dtnitro toluene .  The  identification  of 
the  latter  two  compounds  was  confirmed  by  probe  mass  spectrometry  of 
samples  collected  from  the  HPLC.  A  chromatographic  profile  of  Site  A 
sediment  appears  in  Figure  17  and  the  amounts  of  RDX,  4-A-2,6-B8T,  and  2- 
A- 4, 6 -BUT  appear  in  Table  17.  The  amounts  of  amino-dinf trotoluenes  in  the 
sediment  extracts  were  confirmed  by  GC-ECD  and  GC-MS . 
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Table  17 


CONCENTRATIONS  OF  RDX,  4~A-2,6-DNT, 
2-A-4 , 6-DNT  FOUND  IN  HOLSTON  RIVER  SEDIMENT 


Site  RDX  (ng  g"-1) 

A  358 

B  — 


±±h±^mjLll 

97.3 

18.0 

0.8 

<0.3 


2~A-4,6-DNT(ng  g*1) 

67.5 

3.1 

0.4 


RDX  was  found  only  in  Site  A  sediment  near  the  RDX  lines  1-5  outflow. 
No  partitioning  of  RDX  to  the  sediment  was  observed  further  downstream. 
No  TNT  was  observed  at  any  site  In  the  sediment,  but  the  major  biological 
transformation  products,  the  aminodinltrotoluenes,  were  observed.  The 
amounts  found  Indicate  that  at  least  189  ng  g”1  of  TNT  were  sorbed  at  Site 
A,  24.2  ng  g”1  at  Site  B,  1.4  ng  g”1  at  Site  C,  and  less  than  0.3  ng  g"* 
at  Site  D.  Our  Phase  II  studies  indicated  that  TNT  would  transform 
rapidly  In  the  presence  of  high  microbial  populations  and  added  organic 
nutrient.  Therefore,  it  is  not  surprising  to  observe  this  transformation. 
In  the  aqueous  phase,  however,  we  believe  biologic  transformation  Is  uil 
because  of  the  low  microbial  population. 

I.  Comparison  Of  Actual  RDX  Concentrations  To  Those  Prom  The  Simulation 
of  the  Two-Dimensional  Compartmentalised  River  System 

The  predicted  and  actual  concentrations  of  RDX  found  along  the 
North  Bank,  Center,  and  South  Bank  of  the  Holston  River  appear  in 
Table  18. 
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Table  18 


PREDICTED  AND  ACTUAL  CONCENTRATIONS  (ppb)  OF  RDX  DETERMINED 

IN  THE  HOLSTON  RIVER 


North 

Bank 

Center 

South 

Bank 

Site 

Predicted 

Avg. 

Found 

Predicted 

Avg. 

Found 

Predicted 

Avg. 

Found 

A  (  0.1  km) 

37.0 

55.0 

4.8 

1.0 

0.028 

0.090 

B  (  4.5  km) 

24.0 

20.1 

13.0 

16.0 

3.5 

7.2 

C  (  7.9  ka) 

16.0 

14.1 

12.0 

12.0 

6.0 

6.0 

D  (33.0  ko) 

9.0 

6.5 

9.0 

6.5 

9.0 

6.5 

At  Site  A,  the  found  RDX  concentration  on  the  North  Bank  ia  higher 
than  the  predicted  value  because  a  very  high  SDX  concentration  hugs  the 
North  Bank.  The  higher  predicted  RDX  concentration  in  the  center  of  the 
river  and  the  predicted  lower  concentration  at  the  South  Bank  indicate 
that  the  predicted  lateral  dispersion  is  faster  than  the  actual  rate  at 
Site  A. 

At  Site  B,  the  higher  predicted  North  Bank  RDX  concentration  and  the 
lows*  ‘edicted  Center  and  South  Bank  RDX  concentrations  indicate  that  the 
pre<H  i  lateral  dispersion  is  beginning  to  approximate  the  actual  value. 

However,  at  Site  C,  che  predicted  and  actual  values  agree  quite 
closely,  indicating  that  the  lateral  flow  rate  has  increased  between  Sites 
B  and  C.  This  ie  expected,  since  Site  C  is  much  shallower  than  Site  B 
<1.4  m  versus  2.5  o) . 

At  Site  D,  the  one-dimensional  compartmentalized  river  system  has 
been  employed  for  13  km,  and  complete  mixing  of  the  river  la  predicted. 
The  actual  values  confirm  the  river's  homogeneity;  however,  they  are 
slightly  lower  than  the  predicted  values.  Part  of  the  difference  is 
attributable  to  further  dilution  from  small  input  streams,  surface  runoff 
(rain  occurred  the  night  before  the  sampling  at  Site  D),  and  rain. 


The  actual  values  demonstrate  that  lateral  dilution  is  the 
controlling  factor  In  predicting  the  RDX  concentration  at  locations  near 
the  HAAP.  Once  the  river  reaches  a  steady  state,  we  expect  that  the  RDX 
concentration  will  remain  constant  at  any  time  of  the  day  as  long  as 
environmental  conditions  are  constant.  Dilution,  therefore.  Is  the 
primary  factor  in  controlling  the  RDX  concentration  In  the  Holston  River, 
and  this  phenomenon  can  adequately  he  described  by  the  two-dioens ion&l 
compartmentalized  river  system. 

J.  Comparison  Of  Actual  TOT  Concentrations  To  Those  From  The  Two 

Dimensional  Compartmentalized  River  System 

As  stated  previously,  the  low  discharge  level  of  TNT  and  the  dilution 
of  waetestreams  by  the  Holston  River  prevented  the  collection  of  reliable 
analytical  data  for  model  validation.  However,  since  the  two-dimensional 
compartmentalized  river  system  could  adequately  describe  the  dilution 
processes  in  the  Holston  River,  we  can  compare  a  few  of  the  data  points 
found  at  Site  B  with  those  predicted  by  the  model  simulation. 

The  Table  15  data  indicate  low  (0.05  ppb)  TNT  concentrations  along 
the  North  Bank  (Sites  B-l  to  B-4)  increasing  to  0.17  ppb  near  the  center 
of  the  river.  The  simulation  results  indicate  concentrations  of  1.4  ppb 
on  the  North  Bank  and  0.75  ppb  in  the  center  of  the  river  during  the  time 
of  sampling  (12:00  P.M.).  The  much  lower  found  values  compared  to  the 
predicted  values  on  the  North  Bank  could  arise  from  several  factors.  One 
la  that  the  velocity  of  the  water  on  the  North  Bank  is  much  slower  than 
the  velocity  in  the  center  of  the  river.  On  a  comparative  basis  for 
waters  traveling  between  points  A  and  B,  water  near  the  banks  travels  much 
slower  than  water  in  the  middle  of  the  river.  The  longer  residence  time 
between  the  two  points  allows  more  transformation  to  occur,  thus  producing 
TNT  concentrations  lower  than  those  observed  in  the  center  of  the  river. 
In  the  model  simulation,  average  river  flow  velocity  is  used.  Thus,  for 
accurate  modeling  of  the  North  and  South  Banks,  actual  flow  velocities 
along  these  banks  should  be  used. 


A  second  factor  contributing  to  higher  model-predicted  values  on  the 
North  Bank  Is  the  use  of  an  average  river  depth  value  is  the  calculations 
of  the  photochemical  rate  constant.  Since  the  North  Bank  waters  begin  at 
"zero"  depth  and  increase  to  an  “average"  depth,  within  this  zero-to- 
average  depth  range,  the  photochemical  rate  constant  la  higher  than  the 
value  used  In  the  model  simulation.  Therefore,  the  simulation  model  will 
again  predict  higher  values  than  those  observed.  This  can  be  corrected  in 
the  model  by  using  photochemical  rate  constants  determined  from  shallower 
depths . 

In  the  center  of  the  river,  the  predicted  value  (0.75  ppb)  and  the 
actual  value  (0.17  ppb)  differ  by  approximately  a  factor  of  four.  This 
may  indicate  that  the  model  simulation  is  beginning  to  approximate  the 
transformation  processes  in  the  Holston  River;  more  data  points  are 
necessary  to  confirm  this  expectation.  In  any  case,  a  factor  of  four  may 
be  a  reasonable  value  when  modeling  real  world  environments* 


K.  Persistence  of  SSX,  TAX,  and  HMX  in  the  Holaton  River 

During  the  field  study,  SEX,  TAX,  and  HMX  were  being  discharged 
simultaneously  with  TNT  and  RDX.  The  HMX  concentrations  in  the  discharge 
lines  were  low  (4  to  60  ppb)  and  dilution  from  the  Holston  River  reduced 
the  concentration  of  this  component  to  below  detection  Halts. 

SEX  and  TAX,  however,  were  observed  at  all  sampling  sites.  Their 
concentration  profile  paralleled  that  of  RDX  (Figure  12)  and  indicated 
river  homogeneity  at  Site  D.  From  the  available  data,  it  appears  that  SEX 
and  TAX  will  be  persistent  chemicals  in  the  aquatic  environment. 


V  DISCUSSION 


The  results  of  this  study  indicate  that  the  SRI  model  can  be  used  as 
a  predictive  tool  to  estimate  the  concentrations  of  TNT  and  RBX  in  the 
Holston  River*  The  SRI  computer  model  proved  to  be  reliable  (factor  of 
1.5}  for  estimating  RDX  concentrations  in  the  one-dimensional  compart¬ 
mentalized  river  system  at  distances  beyond  20  km  from  the  discharge  point 
or  at  locations  where  river  homogeneity  has  been  achieved.  However,  due 
to  the  incomplete  mixing  of  the  Holacon  River  at  locations  near  the 
discharge  points,  the  one-dimensional  compartmentalized  river  system  could 
not  be  used  because  a  basic  assumption  within  the  model  considers  the 
river  compartment  a  well-mixed  reaction  vessel.  Since  the  transformation 
processes  that  affect  the  persistence  of  TNT  had  lowered  the  TNT  concen¬ 
tration  below  detection  limits  before  river  homogeneity  could  be  achieved, 
the  one-dimensional  model  was  not  applied  to  TNT. 

Because  of  the  incomplete  mixing  of  the  Holston  River  near  RAAF,  a 
two-dimensional  compartmentalized  river  system  was  developed  that  took 
into  account  the  lateral  dispersion  of  water  across  the  Holston  River. 
When  this  model  was  applied  to  RDX,  the  predicted  RDX  concentrations  in 
the  North,  South,  and  Center  portions  of  the  Holston  River  were  in  good 
agreement  with  the  actual  values.  Since  RDX  is  a  conservative  chemical, 
the  two-dimensional  compartmentalized  river  system  adequately  describes 
the  dilution  of  RDX  in  the  nondlspersed  reach  of  the  Holston  River.  The 
nondispersion  effect  is  dependent  on  the  total  river  flow  velocity  and  the 
SRI  model  can  be  adjusted  easily  to  account  for  this  effect. 

Validation  of  the  SRI  model  using  TNT  was  prevented  by  the  low  input 
concentrations  of  TNI'  during  the  field  study  and  the  poor  dispersion  of 
the  effluent  stream  from  RDX  lines  1-5  into  the  Holston  River.  The  data 
from  the  analysis  of  water  samples  showed  little  penetration  of  TNT  into 
the  middle  of  the  river  and  indicated  that  a  rapid  transformation  of  TNT 
was  occurring  along  the  North  Bank.  The  transformation  rate  appeared  to 
be  more  rapid  than  that  which  could  be  estimated  by  the  model  simulation. 
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This  discrepancy  could  be  accounted  for  by  the  (forth  Bank  waters 
moving  at  a  such  slower  rate  than  the  rate  that  was  measured  at  station  A~ 
1,  25  feet  from  the  shore  (1.78  m3  s"1).  Since  the  two-dimensional  system 
predicted  an  RDX  concentration  of  37  ppb  at  station  A-l  and  212  ppb  was 
found,  we  know  that  the  lateral  dispersion  at  station  A-l  was  greatly 
overestimated.  This  indicates  that  the  RDX  lines  1-5  wastewaters  are 
close  to  the  North  Bank  at  station  A-l  and  are  moving  at  a  flow  between  0 
(shoreline)  and  1.78  m3  a“*.  Thus,  the  residence  time  of  TNT  (and  RDX)  in 
the  North  Bank  compartments  may  actually  be  longer  than  estimated  with  the 
available  data.  Hodel  validation  could  be  achieved  with  TNT  by  allowing  a 
24-hour  continuous  discharge  through  RDX  lines  1-5  and  by  beginning  to 
collect  samples  at  6:00  AM  at  Site  C  followed  by  successive  sampling  at 
locations  downstream.  In  this  way  we  would  obtain  better  analytical  data 
at  higher  concentrations  of  TNT  with  which  to  compare  with  the  predicted 
values . 

We  did,  however,  obtain  sufficient  data  to  estimate  a  transformation 
rate  constant  for  TNT  between  Sites  A  and  B.  RDX  was  used  aa  a  marker 
compound  to  estimate  the  lateral  dilution  factor  between  Sites  A  and  B. 
This  factor  was  calculated  to  be  10,4.  The  initial  TNT  concentration  at 
Site  A  (3.3  ppb)  is  expected  to  decrease  to  0.31  ppb  in  the  absence  of 
transformation*  Since  a  TNT  concentration  of  0.17  ppb  was  found  at  Site  B 
(4.5  km),  we  estimated  that  0.14  ppb  of  TNT  underwent  transformation.  The 
observed  transformation  rate  constant  was  estimated  as  follows: 


lit-— —  -  k  t 

C  p 

lo-~|y  -  kp(0.115  days) 
5.2  day”1-  k 


If  we  assume  an  average  depth  of  1.8  m  between  sites  A  and  B  and  use 
the  thin  film-thick  film  theory  to  correct  for  depth,  the  expected 
photoysis  rate  constant  should  be 
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which  is  in  good  agreement  with  the  measured  value 


These  results  suggest  that  the  laboratory  measurement  of  kp  can  be 
used  reliably  to  estimate  a  value  in  the  river*  The  good  agreement 
between  values  of  kp  also  indicate  that  other  processes  such  as 
biotransforsation  are  not  securing  to  any  extent.  Biotransformation  did 
occur  in  the  sediment;  however,  due  to  the  low  microbial  population  in  the 
river,  this  process  is  not  believed  to  be  occurring  to  any  extent  in  the 
major  portion  of  the  river.  Biotransformation  may  occur  in  isolated 
pools,  but  the  overall  loss  would  not  be  observed  in  the  river. 

That  no  TNT  is  found  in  the  sediments  suggests  that  biotransformation 
is  rapid  in  this  compartment.  How  this  affects  the  sorption  partition 
coefficient  for  TNT  is  not  known  at  this  time.  The  sediment  loading  of 
TUT  transformation  products  does  decrease  between  Sites  A.  and  D,  which 
reflects  the  lower  aqueoua-phase  concentration  of  TNT  reaching  the 
downstream  site. 

A  number  of  assumptions  used  in  this  study  deserve  further  comment. 
First,  we  have  assumed  that  homogeneity  of  the  Holston  River  will  be 
achieved  at  20  km  for  the  given  flow  rate.  The  predicted  and  observed 
values  for  RDX  at  Site  C  are  in  reasonable  agreement,  which  indicates  that 
this  was  a  good  assumption.  However,  the  observed  and  predicted  value  at 
Sites  A  and  B  may  have  been  in  better  agreement  if  lateral  transport  was 
reduced.  Thr.j,  for  predicting  pollutant  concentration,  the  dependence  of 
complete  mixing  (homogeneity)  on  flow  rate  should  be  better  understood  for 
modeling  purposes. 

Second,  we  have  assumed  pseudo-first-order  kinetic  behavior  for  the 
photolysis  of  TNT.  Our  studies  have  shown  that  the  photolysis  depends  on 
other  species  in  natural  waters  (such  as  huaic  acids),  and  the  exact 
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kinetic  dependence  is  unknown,  This  added  complexity  may  cease  some 
variation  in  the  pseudo-first-order  rate  constant  at  different  sampling 
locations.  The  magnitude  of  error  resulting  fro*  this  aaauapti&D  is 
unknown. 

The  enviroameatal  photolysis  rate  constant  was  also  measured  at  soon 
and  was  assumed  to  decrease  to  zero  is  the  sight  and  increase  frost  zero  to 
the  noon  value.  We  assume  that  this  rise  and  fall  of  the  rate  constant 
followed  a  sinusoidal  distribution.  We  believe,  however,  that  deviations 
from  this  assumption  would  produce  an  error  no  greater  than  a  factor  of 
two,  based  os  our  measurement  of  rate  constants  at  noon  and  4:00  ?.M.  in 
our  laboratory. 

We  have  assumed  an  average  depth  for  the  calculation  of  the 
environmental  photolysis  rate  constant.  We  believe  this  is  a  reasonable 
approximation  once  TUT  is  distributed  across  the  river.  However,  in  the 
initial  stage,  the  discharge  is  hugging  the  North  Bank  where  the  depth  is 
much  less  than  it  is  further  Into  the  river.  This  will  tend  to  increase 
the  photo  transformation  rate  constant  for  that  section  of  the  river  until 
average  depth  distribution  is  achieved. 

Xu  conclusion,  this  study  demonstrates  that  models  can  be  used  to 
reliably  predict  BOX  concentration  in  the  aquatic  environment,  and  that, 
therefore,  they  can  be  useful  in  hazard  assessments.  Insufficient  data 
were  obtained  to  make  a  Bimiliar  claim  for  TNT,  due  to  low  input  concen¬ 
tration  during  the  collection  of  field  samples.  However,  from  the  data 
that  were  obtained,  we  conclude  chat  the  simulation  model  will  approximate 
environmental  TOT  concentrations  in  the  Hoi a ton  River.  This  tends  to  make 
the  use  of  pseudo-first-order  rate  expressions  creditable  for  describing 
Che  loss  and  movement  of  aquatic  pollutants  and  suggests  that  the 
laboratory  determination  of  transformation  race  constants  can  be  usefully 
applied  to  aquatic  systems. 

The  study  also  demonstrates  the  much  greater  complexity  of  a  real 
aquatic  systea  compared  to  our  simulated  model  river  and  the  need  to 
incorporate  more  of  the  real-world  complexity  into  the  planning  of  field 
experiments  and  interpretation  of  the  data.  The  great  value  of  having  a 
conservative  chemical  available  with  which  to  trace  out  the  flow 
characteristics  of  the  river  is  evident. 
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